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Abstract. GALFACTS is a large-area spectro-polarimetric survey on the Arecibo 
Radio telescope. It uses the seven-beam focal plane feed array receiver system (ALFA) 
to carry out an imaging survey project of the 12,700 square degrees of sky visible from 
Arecibo at 1 .4 GHz with 8 192 spectral channels over a bandwidth of 300 MHz sampled 
at 1 millisecond. The aggregate data rate is 875 MB/s. GALFACTS observations will 
create full-Stokes image cubes at an angular resolution of 3.5' with a band-averaged 
sensitivity of 90 fiiy, allowing sensitive imaging of polarized radiation and Faraday Ro- 
tation Measure from both diffuse emission and extragalactic sources. GALFACTS is a 
scientific pathfinder to the SKA in the area of cosmic magnetism. Key to magnetism 
science with the SKA is the technique of RM synthesis. The technique of RM synthe- 
sis is introduced and we discuss practical aspects of RM synthesis including efficient 
computational techniques and detection thresholds in the resulting Faraday spectrum. 
We illustrate the use of the technique by presenting the current development of the RM 
synthesis pipeline for GALFACTS and present early results. 



1. Introduction 

Linear polarization of radio sources contains information on magnetic fields in these 
sources, and Faraday rotation of the plane of polarization provides information on the 
direction and magnitude of the magnetic field along the line of sight. As such, ob- 
servations of linear polarization of radio sources provide the most widely applicable 
probe of cosmic magnetic fields on scales from galaxies to clusters of galaxies. Finding 
polarized sources in survey images and fitting their parameters forms the basis of this 
analysis. 

The Galactic ALFA Continuum Survey Consortium is using the Arecibo telescope 
and ALFA to carry out a sensitive, high resolution, spectro-polarimetric s urvey of the 
region of the sky visible with the Arecibo telescope ("Taylor & Salteill2blClf) . A key ob- 



servational objective of the GALFACTS is to image the polarized emission from both 
discrete objects and the dilfuse interstellar medium of our Galaxy and to derive polar- 
ization properties, including Faraday Rotation Measures (RMs) for a vast population of 
extragalactic sourc es. A detailed discussion on the GALFACTS project can be found 
in Taylor & Salter ( lOlOD . 

The objective of GALFACTS is to produce high sensitivity Stokes /, Q, U and V 
multi-frequency images of the polarized emission over the full range of Galactic lati- 
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tudes from the mid-plane to the pole. GALFACTS has a bandwidth of 300 MHz over 
1225 — 1525 MHz, allowing imaging of Faraday rotation measures. The Icr survey 
sensitivity over the full band is ~ 80 fiJy per Stokes parameter {AT„ns ~ 0.8 mK) with 
a 3.4' resolution (maximum baseline 220m); The total sky coverage of the full survey 
is 12,734 degrees which requires a total observing time 1600 hours. Since the GAL- 
FACTS survey produces terabyte sized data sets, the data processing pipeline raises 
considerable challenges. Here we limit our discussion to the RM synthesis methods 
applied to the processed image data, a detailed description of the GALFACTS pipeline 
is given in ,Guram et al, (.2011.) . 



2. Rotation Measure Synthesis 

Polarimetric radio observations enable the study of synchrotron emission radiated by 
relativistic electrons as they are accelerated by magnetic fields. Multifrequency ob- 
servations allow RM Synthesis (Burn| [l966l) to solve for the unknown Faraday depth, 
polarized intensity, and polarization angle simultaneously. The s ource is typically iden- 



tified by the maximum value in the Faiaday spectrum. FoUowing lBrentjens & de Bruynl 
(I2OO5.) the Faraday depth (p is defined as 



</>(r) = 0.81 r -drradm"^, (1) 
Jo 

where B\\ is the line of sight magnetic field component, rig is the thermal electron 
density, dr is an infinitesimal path length, with the integral taken from the observer to 
the point x. 

The complex polarized intensity ^('^^) - Q + iU is the Fourier transform of the 
Faraday dispersion function F(0), 



F{(t))e^'^^^d(p. (2) 



The Faraday rotation measure RM is defined as the slope of a polarization angle 
T versus A^: 

«MU)=-^ (3) 

where 

1 , U 

'F^-tan--. (4) 

Once a polarized source has been detected, the observed polarized intensity p - 
^JQ^ + must be corrected for polarization bias. Since p is a. positive-definite quan- 
tity, the noise in Stokes Q and U results in a positive value for p even if no signal is 
present . The statistics of p with a signal po and noise ctqu is given by the Rice distri- 
bution ( Rice| [l945). This matter is complicated by the uncertainty in the Faraday depth 



of the source. This introduces a stronger bias in the polarized inte nsity than the well 
known polarization bias but can be corrected for (iGeorge et al.ll201ll) . 
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3. Application to GALFACTS 



The first GALFACTS dataset (Nl region) is a cube of 5314 by 1074 pixels with 4096 
spectral channels. The large number of pixels means that to do effective RM synthesis 
is computationally expensive. To solve this RM synthesis can be run in parallel on each 
pixel in the image plane with the results being combined to form one Faraday depth 
cube. This is currently not part of the main GALFACTS pipeline (iGuram et alj 2011 ) 
but w ill be integrated as an application on http : //www . cyberska . org ('Kiddle et all 
201 lb . Another way to approach the RM synthesis is to average the Stokes Q and U 
data in frequency, correctly selecting the window to average so that no information is 
lost in the result Faraday depth cube. For the Nl region RM synthesis was completed 
over a range of -1200 to 1200 rad.m"^ with an averaging of 0.4 MHz in the Q and 
U. To visualize the Faraday depth cube is difficult with such a large dataset and to 
over come this limitation a moment map was created. A preliminary RM map from the 
Nl GALFACTS observations run is shown in Figure 1 with a typical compact source 
Faraday spectrum given in Figure 2. Residual imaging artefacts dominate the compact 
sources but lots of structure can be seen in the diffuse emission. 
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Figure 1: A moment map of a subset of the Nl region showing structure in the 

Faraday spectrum 
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Figure 2: The typical Faraday spectrum of a compact source in the preliminary Nl 

data 
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